Critical Current Oscillations in Strong Ferromagnetic Pi- Junctions 
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We report magnetic and electrical measurements of Nb Josephson junctions with strongly ferromagnetic 
barriers of Co, Ni and NigoFe2o (Py). All these materials show multiple oscillations of critical current with 
barrier thickness implying repeated 0-n phase-transitions in the superconducting order parameter. We show in 
particular that the Co barrier devices can be accurately modelled using existing clean limit theories and so that, 
despite the high exchange energy (309 meV), the large I c Rjv value in the 7r-state means Co barriers are ideally 
suited to the practical development of superconducting 7r-shift devices. 

PACS numbers: 74.50.+r, 74.25.5V, 74.78.Db, 74.25.Ha. 



Although the interplay of superconductivity and ferromag- 
netism has been the subject of study for many decades 1 1], 
theoretical and experimental investigations into the properties 
of superconductor-ferromagnetic metal (S/F) heterostructures 
have seen an upsurge in interest in recent years following the 
experimental observation of to tt transitions in the supercon- 
ducting order parameter in S/F thin films by Ryanzanov et al. 
01 and by Kontos et al. 0|. In terms of the Josephson re- 
lationship I c = Io sin A</>, where Acj) is the phase difference 
between the two S layers, a transition from the to tt states 
implies a change in sign of Io from positive to negative. Phys- 
ically, such a change in sign of Iq is a consequence of a phase 
change in the electron pair wave function induced in the F 
layer by the proximity effect. Experimentally, measurements 
of I c are insensitive to the sign of Io and hence the absolute 
I c is measured; this implies that a change in state from to tt 
will lead to a zero crossing of I c and a sharp cusp at I c = 0. It 
is possible to describe the I c dependence with ferromagnetic 
thickness (dp) by the generic expression 
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where di is the thickness of the ferromagnet corresponding 
to the first minimum and I c i?7v(do) is the first experimental 
value of I c Rn (Rn is the normal state resistance). Transitions 
can be observed as oscillations in the critical temperature (T c ) 
of S/F multilayers |4|, la LZIl as well as oscillations in th e I~ 
of S/F/S junctions with both thickness of the F layer H@H(ill 
and, for weak ferromagnets whose exchange energy (E ex ) is 
comparable to kgT c of the superconductor, with temperature 

The majority of studies of S/F/S structures have concen- 
trated on the use of weak ferromagnets, such as Cu^Nii-a; 
and Pd^Nii-a;, where the temperature dependence can be ob- 
served and where thickness of the ferromagnetic layer over 
which oscillations in critical temperature or current are ob- 
served can be comparatively large. Even where strong ferro- 
magnets have been used, a significant ma gne tic 'dead' layer 
corresponding to a loss in total moment 11211 is usually ob- 
served which complicates the modeling - see table|U 



In the dirty limit where the mean free path L < dp and 
L < hvf/E ex the two decay lengths, £i and £2, take similar 
values and so multiple oscillations of I c are not observed. In 
contrast, in the clean limit where dp < L and L > hvf / E ex 
the decay of the envelope determining the modulation ampli- 
tude (£1) can be much larger than the oscillation period (£2). 
Most previous studies, including those using strong ferromag- 
nets, have been performed in the dirty limit; for practical ap- 
plications, in which large I c Rn values are required in the tt 
state it is vital to develop high quality clean limit S/F/S de- 
vices. A recent report using the ferromagnetic intermetallic 
M3AI shows I c oscillations in the clean limit II 311 . but with 
a very large magnetic dead layer which is not accounted for 
in any phenomenological model and which is likely to make 
practical control of the phase state of the junction difficult. 

Co is a proven device material which can be deposited 
in clean thin film form with accurately controlled thickness; 
however, it has not been previously applied in S/F/S junc- 
tions because the exchange energy was considered to be far 
too large. In this letter, we report for the first time measure- 
ments of junctions containing Co barriers, together with com- 
parative studies of Py and Ni barriers. We show that, unlike 
Py and Ni and the weak ferromagnets, the Co data fits excel- 
lently to clean limit theory. As importantly, the magnetic dead 
layer in the Co is less then 1 nm allowing precise control of 
the phase state of Nb/Co/Nb tt -junctions. 

Nb/Co/Nb, Nb/Py/Nb and Nb/Ni/Nb films were deposited 
on 10 x 5 mm silicon (100) substrates coated with a 250 
nm thermal oxide. Simultaneously to growing 10 x 5 mm 
thin films for patterning, identical 5x5 mm films were de- 
posited for magnetic characterization in a vibrating sample 
magnetometer (VSM). All of the layers were deposited by 
d.c. magnetron sputtering at 1.5 Pa and the deposition sys- 
tem was baked out for seven hours and subsequently cooled 
with liquid nitrogen for at least one hour prior to the depo- 
sition, which gave a base pressure better than 5 x 10 -6 Pa 
and an oxygen partial pressure of less than 3 x 10~ 9 Pa. 
The deposition rates are: 1.2nmmin _1 for Co, l.&nmmivr 1 
for Py, QAnmmin~ l for Ni, 2Anmmin~ 1 for Cu and 
YlfSrimminT 1 for Nb. In a single deposition run, multi- 
ple silicon substrates were placed on a rotating holder which 
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passed in turn under three magnetrons. The thickness of each 
layer was controlled by setting the angular speed at which 
the substrates moved under the respective targets and by set- 
ting the target power. When depositing Co, Py, and Ni bar- 
riers, an acceleration curve was programmed which allowed 
the angular speed of the substrates to change monotonically 
as the substrates passed under the relevant targets. The thick- 
ness of the Co (dco), Py (dp y ) and Ni (d/vi) was then depen- 
dent on the substrate position, 0, on the rotating holder. With 
knowledge of the deposition parameters, the rotation was pro- 
grammed such that d(dp)/dd was a constant. This method 
of varying dp guaranteed in all cases that the interfaces be- 
tween each layer prepared under the same conditions and the 
only variation was the thickness. To confirm control over the 
thicknesses we performed X-ray reflectivity of Nb/Co/Nb thin 
films where the Nb layers had a thickness of 5 nm and the Co 
barrier thickness was varied from 0.5 nm to 5.0 nm. A series 
of low angle X-ray scans were made and the thickness of the 
Co layer (dco(observed)) extracted by fitting the period of the 
Kiessig fringes using a simulation package. It was found that 
our expected thicknesses, dc ( expected) > was wel1 correlated 
with d Co ( obaerved } with a mean deviation of 0.2 nm. 

To assist in locating the barrier layer in subsequent focused 
ion beam (FIB) processing, a thin (20 nm) normal metal layer 
of Cu was embedded in the 250 nm thick Nb electrodes lo- 
cated 50 nm from the ferromagnetic layer (a distance greater 
then the coherence length of Nb). At 20 nm, the Cu layer is 
much thinner than the coherence length in Cu and is there- 
fore totally proximitised by the Nb and plays no part in the 
electrical properties of the junctions. The thicknesses of the 
F barriers for the Co and Py junctions was varied from ap- 
proximately 0.5 nm to 5 nm and the Ni barrier thickness was 
varied from 1.0 nm to 10 nm. The films were patterned using 
optical lithography, followed by broad beam Ar ion milling 
(3 mAcm -2 , 500 V beam) to produce micron scale tracks 
and contact pads, to allow four-point measurements to be per- 
formed. These tracks were then patterned with a Ga + FIB 
(Philips/FEI FIB 200) to achieve vertical transport with a de- 
vice area in the 0.2 — 1 /im 2 range. The FIB technique for pro- 
cessing such junctions is described in detail elsewhere 11 411 . A 
micrograph of one of the junctions processed using the FIB for 
this work is shown as the inset of Fig. EJa). 

To investigate the magnetic properties of our films we have 
measured, using a VSM at room temperature, the magnetic 
moment per unit surface area of the films as a function of dp 
(see Fig. 0. Extrapolating the least-squares fit of the data 
in Fig. Ogives magnetic dead layers of ~ 0.75 nm for Co, 
~ 0.5 nm for Py and ~ 1.5 nm for Ni. The causes of a mag- 
netic dead layer are attributed to factors such as lattice mis- 
match causing elastic deformation, formation of amorphous 
interfaces and a breakdown in the crystal structure leading to 
a reduced exchange interaction between neighbouring atoms 
and hence a reduction in Tcurie and E ex I15l lla Il7ll . In the 
case of Co and Ni, both the thickness of the dead layer and the 
total moments for a given thickness greater than the dead layer 
are close to those reported in systematic studies Nb/F bilayers 
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FIG. 1: (Color online) Scaling of the magnetic moment per unit area 
vs. Co, Py and Ni barrier thickness at room temperature. Inset: hys- 
teresis loops for Co and Ni at room temperature. 




Transport measurements were made in a liquid He dip 
probe. The differential resistance as a function of bias current 
of the junctions was measured with a standard lock-in tech- 
nique. The I c was found from the differential resistance as 
the point where the differential resistance increases above the 
value for zero bias current. Rm was measured using a quasi- 
dc bias current of 3-5 mA, this enabled the nonlinear portion 
of the I-V curves to be neglected, but was not large enough to 
drive the Nb electrodes into a normal state. A dV/dI(I) and 
V(I) plot for a typical Nb-Py-Nb Josephson junction is shown 
inset of Fig. |2jb). In general, the critical currents of the de- 
vices measured for this work ranged from 500 /iA to below 
the minimum sensitivity of our apparatus (50 nA), while the 
normal state resistances were in the range 1.0 mfl to 100 mf2. 

The variation of I c Rn as a function of Co, Py and Ni thick- 
ness at 4.2 K is shown in Figs. |2ja-c) respectively. Each point 
in these Figs, corresponds to the mean of several junctions 
with different areas; the vertical error bars are derived from 
the measured variation in I c Rn and a small noise contribu- 
tion due to the current source. From X-ray reflectivity results, 
as discussed above, we take the error in dp for all F barriers 
to be w 0.2 nm. All of the devices shown in this data set 
presented Shapiro steps upon the application of microwaves 
and an I c modulation with applied field Ha- For the devices 
where dco > 5 nm, dp y > 5 nm, and d^i > 13 nm some re- 
duction of the differential resistance around 7 = was seen, 
but did not show a measurable supercurrent at 4.2 K. As ex- 
pected, I c Rn for the Co, Py and Ni decreases exponentially 
and in an oscillatory fashion with dp. 

The Co data was modeled using Eq. Q; from the model fit 
shown in Fig. |2ja) we find that the period of the Co oscilla- 
tions is ~ 1.9 nm, hence £2 — 0.3 nm and £1 ~ 3.0 nm. This 
gives a £2/^1 ratio of ~ 0.11. A full theoretical treatment in- 
volves solving linear Eilenberger equations 1 19] and gives Eq. 

tanh - A. = % (2 ) 
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FIG. 2: (Color online) Characteristic voltage as a function of Co, 
Py and Ni thickness at 4.2K. The solid lines are fits to Eq. 0, the 
dotted lines are fits to Eq. J5j, and the dashed line is a fit to Eq. (3 
as described in text. Inset (a): An FIB micrograph of a typical Nb- 
Co-Nb Josephson junction. Inset (b): V(I) and dV/dI(I) plotted 
for a Nb/Py/Nb Josephson junction at 4.2K. Inset (c): normalised 
magnetization as a function of warming and cooling temperature for 
a Nb/Ni/Nb trilayer where djvi — 9 nm. 



where £ e ff is the effective decay length given by £ e y^ = 

1 + i£j , £ is the Ginzburg-Landau coherence length 
and £jj is a complex coherence length. In the clean limit 
l + L^- 1 » imax{ln(l + Lf^lnfLQ 1 )}. The solu- 



tion of Eq. (0 gives 
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and the numerical solution is shown in Fig. (0. 

Following the method of Gusakova and Kupriyanov lll9ll 
we find from Fig. (|3j that the experimental ratio £,2/^1 — 
0.1 corresponds to two inverse magnetic lengths of L/^jj — 
16.5 and L/£jj ~ 18.7. By assuming L/£ ~ 0.1 and for 
the estimated parameters £1 ~ 3 nm and £ 2 — 0.3 nm we 
obtained, from the inset in Fig. that for L/£jj ~ 16.5 




FIG. 3: (Color online) The dependence of £a/£i with inverse mag- 
netic length, L/£h, calculated for different ratios of L/fo. Inset: 
inverse decay length, ~ /(L/fo) for when L/£h — 0.1 



and L/£h — 18.7, L ~ 5 nm. Inputting these values into 
Eq. gives vf{C6) ~ 2.8 x 10 5 ms" 1 which is similar to 
other reported values of vf(Co) J20I1 and E ex ~ 309 meV. 
To validate the mean free path of our Co thin film we have 
estimated L(Co) in a 50 nm thick Co film by measuring its 
resistivity as a function of temperature using the Van der Pauw 
technique \ 2l}. Following the method described by Gurney et 
al. we estimate that our Co has a mean free path of L(C6) ~ 
10 nm which justifies our assumption of a clean limit fit. As a 
comparison we have also modeled the Co oscillations with a 
simpler theoretical model given by Eq. 1 23] 



I C R N oc 
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2E ex dF/hvf 

As in the case of Eq. Q the fitting between the theoretical 
model and the experimental data is good (see Fig. |2|(a) dashed 
line) and in particular the best fitting has been obtained by 
using v F = 2.8 x 10 D ms -1 and E ex = 309 meV. 

In contrast, the Py and Ni data cannot be fitted entirely in 
the clean limit. In the case of Py, Eq. @, closely matches the 
experimental data up to a thickness of ~ 3.6 nm and in the 
case of Ni the oscillations follow the clean limit theory to ~ 7 
nm. Above these values a better fit is obtained using a formula 
for a diffusive and high E ex F 1 24]: 
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where A is the superconducting energy gap, oj m is the Mat- 
subara frequency and is given by uj m = 7rTfcs(2r7i+l) where 
T is the transmission coefficient and m is an integer number. 
k u = (1 + 2 I oj m I t/Ti) — 2iE ex T/h and fi — cos 8 where 
8 is the angle the momentum vector makes relative to the dis- 
tance normal to the SF interface. L is given by VfT and r is 
the momentum relaxation time. For Eq. @ the only fitting 
factor, besides the numerical prefactor, was the strength of the 
exchange interaction (E ex /hvf). In the case of Eq. (|3J a suit- 
able Vf, A and E ex had to be chosen. To fit Eq. to Py and 
Ni data we used: Vf(Py) — 2.2 x 10 5 m/s and L Py ~ 2.3 nm, 
and Vf(Ni) = 2.8 x 10 5 m/s and L Nl ~ 7 nm and A = 1.3 
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Ref 


1.2 


1.6 


Ni 2 oFe 8 o 


2.2 xlO 5 


95 


0.7 




1.4 


0.46 


Ni 2 oFe 8 o 


2.2 xlO 5 


201 


0.5 
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1.8 


2.0 


PdgoNiio 


2.0 xlO 5 


35 


_ 
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1.7 


1.0 


Ni 


2.8 xlO 5 


200 




UPj 


2.3 


0.86 


Ni 


2.8 xlO 5 


107 




[81 


4.1 


1.2 


Ni 


2.8 xlO 5 


80 


1.5 


* 


4.6 


0.45 


Ni 3 Al 


1.5 xlO 5 


86 


5-8 




3.0 


0.3 


Co 


2.8 xlO 5 


309 


0.8 


* 



TABLE I: A summary of reported parameters for different material 
systems. DL stands for 'Dead Layer Thickness' and * 'This paper'. 



meV. These values are consistent with the ones used in Eq. 
@ and elsewhere JHU. From the oscillations period we can 
estimate the E ex of the Py and Ni barriers - from Eq. ® the 
periodicity is given by L osc ~ hvf/2E ex and hence we can 
deduce E ex . The exchange energies for Py and Ni were found 
to be 201 meV and 80 meV respectively. E ex (Py) is approxi- 
mately double that measured in Nb/Py/Nb junctions deposited 
with epitaxial barriers where E ex was found to be 95 meV @] • 
E ex (Ni) is close to other reported values by photoemission 
experiments [2o|. The smaller than expected E ex is a conse- 
quence of impurities and possibly interdiffusion of Ni into Nb. 
In the case of Ni, we have measured the magnetisation satu- 
ration as a function of temperature (M (T)) (see inset of Fig. 
|2c)) so that an estimate of the Ni Curie temperature could 
be made. This provides further evidence that our Ni data is 
consistent and that the Ni isn't grossly degraded. To be sure 
our data was not influenced by thermal diffusion during this 
measurement we measured M(T) for both the heating and 
cooling phases and it is shown that the cooling curve follows 
the warming one. The warming data is modeled by the for- 
mula: M(T)/M(0) = (1 - T/T clme y j where M(0) is the 
magnetisation at absolute K, T is the measuring tempera- 
ture, and (3 and Tcurie are fitting parameters. We estimate 
T cur i e ~ 571 K which is in agreement with Curie tempera- 
ture measurements of Ni in S/F bilayers [25]. 

The fully clean limit behavior observed with the Co barriers 
arises most obviously from the use of a pure element, but also 
from the vertical coherence likely Il27ll even in noncrystalline 
heterostructures. The high E ex results in a short oscillation 
period implying a need for A-level control of layer thickness 
for practical devices; however Co and Co-alloys form the ba- 
sis of current spintronic device production, and precision layer 
control and excellent compatibility with tunnel barriers 12811 
have been demonstrated in many industrial processes 112911 . Co 
is an attractive material for qubits and other novel devices. 

In summary, we have measured critical current oscillations 
in Co junctions as a function of Co barrier thickness which 
indicates that the devices are strongly in the clean limit. This 
results in a higher I c Rn values in the it state compared to 
I c Rn values in the dirty limit. We also present complemen- 
tary critical current oscillations through Py and Ni barriers. 



The oscillations in I c Rn with dp are indicative of — ir 
crossovers, but without conducting phase-sensitive measure- 
ments the specific phases of the individual junctions cannot 
be known. The oscillations also showed an excellent fit to 
theoretical models. We have also estimated, from the peri- 
odicity of the oscillations, the exchange energies of the Co, 
Py and Ni barriers to be 309 meV, 201 meV and 80 meV re- 
spectively. Results within this paper are summarized in tabled 
alongside results reported elsewhere. Our results are not only 
interesting in their own right, but are a vital experimental step 
towards understanding the physics of quantum electronic de- 
vices based on superconductors and are of considerable value 
to the development of quantum information processing. Our 
devices are precursors to practical implementations into qubits 
and other applications in controllable and scalable supercon- 
ducting quantum electronic devices. 
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